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ABSTRACT 


Shock  relations  are  established  for  the  dense  high  temperature 

I  0^ ’’K  ^  *  Introducing 

relativistic  concepts.  It  Is  assumed  that  the  self-consistent  static 
electric  field  does  not  play  a  significant  role  In  the  shock  relations. 
From  these  relations,  the  heating  effect  due  to  the  passage  of 
the  shock  waves  through  a  degenerate  gas  Is  deduced  explicitly.  In¬ 
dicating  large  heating  effects.  The  results  may  be  of  significance 
In  understanding  the  supernova  phenomenon. 
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I.  INTRODUCTION 


The  Rankine-Hugoniot  relations  for  shock  waves  are  the  conservation 

relations  for  mass,  momentum  and  energy  of  the  flow  in  a  continuous  medium* 

These  have  long  been  known  for  the  case  of  the  ideal  polytropic  gas.^^^  P.  Lai 

(2) 

and  P.  L.  Bhatnagar  have  indicated  one  method  for  extending  these  relations 
to  the  degenerate  Fermi-Dirac  gas.  In  this  paper,  we  shall  rederive  these 
relations  for  the  F-D  gasj  we  find  that  our  results  differ  from  those  of  Lai 
and  Bhatnagar,  because  the  basic  physical  picture  on  which  the  derivation  is 
based  is  itself  different. 

Due  to  the  extremely  high  mass  density  in  the  interior  of  a  white 
dwarf  or  a  presupernova  star  (  ~  10^ 10^  gr./c«c.),  the  matter  is  com¬ 

pletely  ionized  and  highly  degenerate  even  at  temperatures  of  the  order  of 
10^  °  K.  The  electron  gas  plays  a  dominant  role,  and  it  has  been  shown  by 
Landau^^^  that  an  ordinary  sound  wave  cannot  propagate  in  a  neutral  degenera¬ 
tive  medium  unless  the  tenperature  (in  energy  units)  is  cottparable  to  the 
Fermi  energy.  On  the  other  hand,  Landau  has  shown  that  in  the  presence  of 
interaction,  a  new  propagation  mode,  termed  zero  sound,  can  exist. 

It  is  our  purpose  to  establish  the  shock  relations  for  the  dense 
degenerate  gas  and  to  extract  some  of  the  physical  features  a  shock  wave  must 
have  in  such  a  plasma.  We  do  not  consider  here  how  such  a  shock  wave  may  be 
formed,  but  assuming  that  shock  waves  of  the  zero-sound  type  exist,  we  inves¬ 
tigate  their  properties. 

We  make  the  following  approximation.  The  self-consistent  static 
electric  field  is  not  included  in  the  distribution  function  explicitly.  >fe 
expect  that  the  electrostatic  energies  in  the  front  and  rear  regions  of  the 
shock  wave  are  of  the  same  order  of  magnitude  so  that  the  contributions  to 
the  enthalpy  due  to  interaction  energies  will  tend  to  cancel  out.  This 
point  should  be  studied  and  clarified  in  the  future. 
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The  phase  velocity  of  the  zero-sound  is  of  the  order  of  the  velocity 
of  light,  for  the  temperatures  and  densities  considered,  and  this  property, 
together  with  the  conservation  of  the  energy,  gives  a  large  heating  effect, 
perhaps  as  large  as  a  factor  of  20  in  temperature »  This  heating  process  can 
occur  almost  instantaneously  due  to  the  large  propagation  velocity  of  the  zero- 
sound. 

Since  the  existence  of  zero-sound  depends  on  an  interparticle  interac¬ 
tion,  we  have  examined  the  question  of  whether  such  an  interaction  alters  the 
form  of  the  Rankine-Hugoniot  relations.  This  study  is  presented  in  the  Appen¬ 
dix,  where  it  is  concluded  that  the  form  of  these  relations  is  unchanged. 
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II. 


We  shall  assume  extreme  conditions  of  density  and  temperature,  cor¬ 
responding  to  a  presupernova  core..  We  assume 


9  / 

J  ^  10  gr./c«Co 
T  -  10^  ®  K. 


density  of  mass 
Temperature 


(1) 


Under  these  circtimstances,  the  gas  is  completely  ionized,  but  the  electron  gas 
is  strongly  degenerate.  The  ion  distribution  is  principally  determined  by  the 
requirement  of  electrical  neutrality. 

The  conditions  in  (l)  can  be  expressed  most  conveniently  in  terms  of 
new  variables s 


A  ( f.  “  Fermi  momentum) 

'ThC  “ 


(  ^  «  Boltzmann  constant) 


and 


E  2  J<3L 

/K,  -  (  Ej.  =  Fermi  energy) . 


These  parameters  have  the  following  numerical  values: 


(2) 


The  symbols  m  and  c  have  their  usual  meanings;  mc^  is  the  rest  energy  of  the 
electron. 

Following  S.  Chandrasekhar^^^,  one  finds  for  the  free  electron  gas 

the  pressure  P,  number  N  ,  and  internal  energy  U  s 

6  6 


7.  Sex')  -* 


(3) 

a) 

(5) 


with 


=  iL 

^li)=  ?  if-l]'  f(pe i , 


V  refers  to  the  total  ▼olume.  f(i)  can  be  approximated  as 

f(x)  -  2x^  -  3x^  for  X  >>  1 


(6) 

(7) 

(8) 


Compared  to  the  above  values  for  the  electron  gas,  it  can  be  shown 
that' pressure  and  internal  energy  of  the  ion  gas  are  negligible.  Furthermore, 
the  number  of  electrons  is  related  to  the  density  of  the  gas  as 

■  fr  ' 

where  ^  is  mean  molecular  weight  per  electron,  and  H  is  mass  of  the  hydrogen 
atom.  The  mass  density  is  given  by 


S=^H 


Sv-m. 


•i-  IT  t 


J 


(10) 

Having  in  mind  that  x  «  8  for  the  situation  under  consideration,  one  finds  that 
the  internal  energy  of  the  gas  per  unit  mass  can  be  expressed  as 


■M'-f  fW  + 

^X)  t  •fTr'-i'  *(’('=»  O* 

P 

It  follows  that  the  enthalpy  of  the  gas  of  unit  mass,  is  U  ♦  -  ,  can  be 
expressed  as 

.  __  ^xXix>\T- Q 

^  S  ^ 


(11) 


(12) 


Regarding  the  gas  as  polytropic  with  polytropic  index  JT  ,  one  finds 
(putting  €, »  0) 

y  X^-  y  X  ^  Z  »  X  (13) 

*  ''  6x*-  -f'X  7  ^ 

In  the  limit  x  — ►  oo  ,  one  finds  Y  approaches  l*/3.  It  is  known  that  in  the 
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limit  X — ^  0,  for  which  (13)  is  not  valid,  ^  approaches  ^3*  Therefore,  one 
may  state  that  the  degenerate  electron  gas  behaves  approximately  as  a  polytropic 
gas  for  ^ich  }f  varies  between  ij/3  and  5/3,  depending  on  the  density,  and 
weakly  on  temperature. 

The  entropy  of  the  gas  is  given  as 

S  -  i  (U  +  PV  -  NF),  (li4) 


with  F  denoting  the  chemical  potential.  From  the  definition  of  the  chemical 
potential 


F 


dN 


» 


(15) 


we  deduce  that 

M,  P.-HF,  .  -NIL,  (16) 

ion  ion  e  electron  e  r  *  '  ' 

because  ion  and  electron  gas  are  bound  to  each  other  to  satisfy  the  requirement 
of  electrical  neutrality.  Vs  also  find 

\on  «  ®el  ' 


so  that  the  entropy  of  the  gas  can  be  represented  simply  by  the  entropy  of  the 
electron  gas.  This  is  found  to  be  (still  neglecting  £  ) 


or 


— TP - , 

c-»  n'lcr  Je 

me'-  -X^ 


(per  unit  mass) 


(17) 


One  normally  defines  the  sound  velocity  as 


sound 


(|L) 

entropy 


(18) 


Prom  eq.  (17)»  the  entropy  of  the  gas  is  a  function  only  of 


Thus, 
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if  S  is  held  constant,  it  follows  that 

-  0, 

L  >?  /  entropy 


(19) 


since  P  is  a  function  of  x  onlys 


P 
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(20) 


Thus,  we  derive  the  fact  that  ordinary  sound  will  not  propagate  through  a 
degenerate  gas  at  zero  temperature.  This  situation  was  pointed  out,  by  Landau, 
in  discussing  the  zero-temperature  Fermi  liquid.  In  the  presence  of  interac¬ 
tion,  Landau  proved  that  there  exists  another  mode,  called  zero-sound,  which 
propagates  essentially  with  Fermi  velocity. 
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III. 

Let  us  new  consider  a  one-dimensional  shock. 


± 

W/ 

X 

% 

Subscripts  1  and  2  refer  to  the  front  and  rear  regions  of  the  shock,  and  q 
denotes  velocity  of  the  gas  relative  to  the  shock  region.  From  Mf^ler's^^^ 
work,  one  can  write  the  following  equations  expressing  the  conservation  of 
number  of  the  particles  in  the  flov>  aomentum  and  energy  in  the  front  and 
rear  regions  of  the  shock  wave  with  respect  to  the  coordinate  system  to  which 
the  shock  region  is  at  rest! 


t  0^« 

u,  +  A 


*  c*’ 

r».  c'' 


(partielo  number) 


(momentum) 


(energy) 


tnr 

Equation  (23)  can  be  rewritten  as  follows 


(21) 


(22) 


(23) 


(23a) 


These  relations,  expressing  basic  physical  conservation  laws,  must  of  course  be 
satisfied  by  the  degenerate  gas  as  well  as  the  ordinary  gas.  The  relativistic 
conception  is  introduced  in  order  for  this  formalism  to  be  consistent  with  the 
fact  that  the  phase  velocity  of  the  zero-sound  type  shock  wave  is  of  the  order 
of  the  velocity  of  light. 
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The  enthalpy  of  the  gas  per  unit  mass,  i,  in  eq.  (12)  may  be 
rigorously  to  the  lowest  order  in  £,  ,  as 


-h  Tf^'x 


with 


or 


'7>16  * 


MH 


(x-i  -tw’k.  ) 


From  eq>  (2l),  (22),  and  (25),  one  finds  that 


^  xri  ^  ^ 

The  right  hand  side  can  he  approximated,  on  the  aeaunptlon  that  3^, 
values  given  hy  eq*  (2),  as 


'TU, 


OUk 

Mtf 


i(TCy-x,)f 


tb  then  obtain,  solving  for  the  temperatures 


So  long  as  3C2"’^  ^  obtains 


3  CXw-x,; 


It  will  be  shown  later  that  Xg-x^  ^  0  is  valid  always  for  the  physical 
cess.  This  leads  to  the  conclusion  that 


written 


(2U) 


(25) 


(26) 

have  the 


(27) 


(28) 


(29) 

pro- 
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1‘f  T,  ,  (30) 

with  Xg-x^  ^  implies  that  shock  wave  heats  the  degenerate  gas  and 

raises  its  temperature  by  a  factor  of  10  or  more  for  the  case  with 

1«02  It  is  interesting  to  note  that  a  slight  fractional  increase  in 
the  nunler  density  can  produce  a  large  teirperature  rise  in  the  gas. 

Ito  now  show  that 

<V*i>  >  ® 

for  the  physical  process*  from  eq.  (g2)  and  (23a),  one  obtains 


Making  use  of  eq.  (3)  and  (10),  eq*  (31)  can  be  rewritten  as 


or,  from  eq*  (22), 

Solving  for  q  in  eq.  (21),  one  has 


2 

Substituting  q£  in  eq*  (3U)  into  eq.  (33),  one  finds  that 

(  ix.-  -I-  + 

_ _ 4C%. 


Assuming  that 


r 


Ui-hC^)  C3.X,) 


i 


1 


* 


(31) 

(32) 

(33) 

(3U) 

(35) 


the  right  hand  side  of  eq.  (35)  can  be  simplified  as 
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or,  by  eq.  (27), 


2%,  'n,  L  ^ 

Thus,  eq.  (35)  becomes  as  follows: 


^>1  ;««.  21, 


^  21i  _l 

7l>^  a,yx  'Vt 

Consider  first  the  case  where 


(36) 


i)  *2  >  *1  * 

Equation  (36)  appears  to  be  true  because  both,  sides  of  the  equation  become 
larger  than  one  and  positive. 

For  the  case  where 

ii)  X2<  Xj^  , 

one  finds  that  the  left  hand  side  is  positive  and  less  than  one,  while  the 
right  hand  side  of  eq.  (36)  becomes  negative,  yielding  a  contradiction.  More 
careful  study  reveals  that  when  n^  -  Ug,  eq.  (36)  may  hold  true.  However, 

^1  *  ^2  where 

\  ■  % 


Tdiich  is,  of  course,  unphysical,  so  that  it  should  be  discarded.  Thus,  one  is 
led  to  the  condition 

Xg—  ^  0  for  the  physical  processes. 

These  results  imply  that  the  passage  of  the  shock  through  the  medium 
increases  the  density,  and  greatly  increases  the  temperature.  The  implications 
for  advanced-stage  thermonuclear  reactions  seem  clear. 
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IV.  CONCLUSION 

It  is  found  that  the  quantum  nature  of  the  neutral  Fermi  gas  must 
be  taken  into  account  in  considering  the  propkgation  of  an  acoustic  wave 
through  the  gas.  Ve  find  that  a  large  heating  effect  is  to  be  e^ected^ 
which  agrees  in  principle  with  current  ideas  relating  to  supernova  explosions. 
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APPENDIX 


The  Rankine-Hugoniot  relations  are  known  to  be  jump  conditions  for 
shocks,  v^iich  can  be  derived  by  applying  the  three  general  principles  — 
conservation  of  mass,  of  momentum,  and  of  energy  —  to  a  column  of  gas  in  a 
tube.  It  seems  appropriate  to  begin  with  Boltzmann's  equation  as  a  basic 
equation  to  describe  the  transport  phenomena  of  degenerate  gas. 

Let  us  assume  the  column  covers  at  time  t  the  interval 
^^(t)  ^  X  ^  where  a^(t)  and  a^^Ct)  denote  the  positions  of  the  moving 

particles  that  form  the  ends  of  the  column,  and  the  flow  is  supposed  to  be 
continuous  at  the  ends  of  the  column.  Furthermore,  we  assume  that  in  the 
moving  column  there  is  a  point  of  discontinuity  whose  coordinate  x=  ^  (t) 
moves  with  velocity  ^(t)  -  u(t). 

Following  Kadanoff  and  Baym'  ',  one  has  the  following  conservation 
laws  for  degenerate  gas: 


where  ^n(x,t)^  is  the  number  density  of  particles  or  quasi -particles* 
averaged  over  the  various  momentum}  that  is 

and  ^j(x,t)^  is  the  current  density  of  particles  averaged  over  momenta. 


(5) 


« 


Ms  will  use  the  terms  particle  or  quasi-particle  interchangeably,  depending 
on  Tdiether  the  medium  is  a  Fermi  gas  or  a  Fermi  fluid. 


with  E(p,x,t)  total  energy  of  a  quasi  particle.  Here  m  is  a  mass  of  the 
particle,  U  is  the  external  potential,  and  J(x,t)  is  the  momentum  stress 

■SCSj 

tensor,  ejqsressible  as  follows: 


With  £.(x,t),  the  energy  density,  given  by 

^  ^ 

The  fvinction  f(p,p^jx,t)  is  defined  by 


(6) 


(7) 


ftf,  p';  t,t:)  =  3  (g) 

S-n.LpiK.ti'' 

which  can  be  interpreted  as  the  increment  of  particle  energy  at  (p,x,t)  due  to 
the  addition  of  a  particle  at  (p^,x,t).  Finally 

^ u,t:>  s  B Yp  f  (hx, ir)  -n  i)  (9) 


Considering  the  limiting  process  where  the  length  of  column  approaches 
zero,  one  has  generally  equations  of  the  form 

J  \ (if  V*y  )  -t  =0 

A,  ^  ^  '  '  J 

As  long  as  (x,t)  does  not  behave  like  a  delta-function,  one  has 


il 

while  the  second  term  in  eq.  (lO)  yields  always 

•  yl '  =  fcA,)-  fcA,) 

\  /w  , 

The  third  term  in  eq.  (lO)  yields 

/.  r^' 

* 


(11) 


(12) 
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( 


provided  also  that  does  not  behave  like  a  delta  function.  Therefore, 
with  the  assumption 


l£l^  M’?!  »  ,  -for 

eq.  (10)  gives 

fca,}  =  yot^} 


(13) 


(li») 


This  is  the  jump  condition  we  seek. 

With  the  above  reasoning,  one  gets  from  eqs.  (l),  (2)  and  (3),  the 
following  relations,  respectively: 


(1?) 

<Tlj(Vt)>  ■ 

(IT) 

Here  a^(t),  a^(t-)  refer  to  the  points  front  and  rear  of  the  discontinuous 
shock  region.  We  may  also  derive  from  eq.  (l),  together  with  the  second  law 
of  thermodynamics, 


(a^,t)'^  ^  (18) 

where  s  is  the  entropy  per  particle.  Consequently,  in  order  to  obtain  the 
Rankine-Hugoniot  relations  for  a  degenerate  gas,  one  has  to  evaluate  ^  j(x,t)^  , 
^Tj^j(x,t)^,  ^j^  (x,t)^  and  ^ j (x,t)s(x,t)^  ^  explicitly  under  the  specific 
conditions  and  interactions  of  the  gas. 

In  the  case  of  local  equilibrium,  one  has 

j  -  v(x»t)  n(x,t)  ,  (19) 

T^j  -  inVj^(X,t)  Vj(^,t)  n  +  P(^,t)  (20) 


and 
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j  -  vJ|  n  +  (21) 

^  ^  ) 

Tjhich  ol?viously  reduce  to  the  conventional  form  of  the  Rarikine-Hugoniot  relations. 

Thus,  one  is  led  to  the  conclusion  that  the  Rankine-Hugoniot  relations 
have  the  same  structure  for  a  degenerate  gas  as  for  classical  continuous  media, 
but  that  the  terms  should  be  e^qpressed  properly  for  the  degenerate  gas  using 
the  specific  interaction  between  the  particles. 


17 


I  *  '  References 

Courant-Friedrichs .  “Supersonic  Flow  and  Shook  Waves,"  pp.  121f. 

Interscience  puh*  (19U8)* 

L.  P*  Kadanoff  and  G.  Baytn.  "Quantum  Statistical  Mechanics,"  Chapt.  6,  10,  11» 
Benjamin,  Inc.  (1962). 

L.  D.  Landau.  J.E.T.P.  30,  1058-1061*.  (1956). 

L.  D.  Landau.  J.E.T.P,  32,  59-66.  (1957). 


